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We propose a simple alternative proof of the Wronskian representation formula obtained by Gail- 
lard and Matveev for the Darbonx-Poschl- Teller (DPT) potentials. It rests on the use of singular 
Darboux-Backlund transformations applied to the free particle system combined to the shape in- 
variance properties of the DPT. 



PACS numbers; 



I. INTRODUCTION 



In 2002, Gaillard and Matveev established an interesting connection between the Darboux-Crum dressing formula 
and the Darboux-Poschl- Teller potentials. More precisely, they proved the following theorem: 

Gaillard- Matveev's theorem: The trigonometric Darbonx-Poschl- Teller potential (TDPT) with integer parameters 
m > n defined on x G ]0, 7r/2[ by 

, m(m -\- 1) n(n -i- 1) , . 

Vix;m,n) = + ^ \ (1) 

sm X cos^ X 

admits the following Wronskian representation 

Vix;m,n) ^ -2{\nW {ui,...,urn \ x))" , (2) 
where W {ui, Um \ x) denotes the Wronskian of the functions ui (x) , Um {x) which are given by 

Uk {x) = sin(/a;), if 1 < / < rn — n , , 

Um-n+i {x) = sin((TO — n + 21) x), if 1 < / < n. 

Gaillard and Matveev underlined that this result traduces the existence of a chain of Darboux transformations 
connecting the constant potential to the TDPT V{x; m, n). If such a construction was already envisaged in Darboux's 
seminal worksi, the theorem above precises what is the correct choice of "seeds" eigenfunctions of the free Hamiltonian 
. !^ ! to reach this last potential. 

. The proof given in^ and'^ rests on direct evaluations of the Wronskian W (ui, Um \ x) and is quite sophisticated 
5^ ' and technically involved. In this letter, we show that this result can be obtained in a shorter way directly by following 
^ , Darboux's original idea. In a modern formulation the proof makes use of shape invariance and SUSY QM arguments 
associated to singular Darboux-Backlund Transformations (DBT) built from excited eigenstates. If the use of singular 
DBT has already been envisaged in several papers^"—, they are generally considered of less interest, although recently 
proven to be crucial in the study of some two-dimensional superintegrable systems^. 

We also examine the case of the Bessel potentials, which is the first example envisaged by Darboux^. This is a 
confluent case which necessitate to employ Matveev's generalized Wronskian leading then to a Wronskian version of 
the Rayleigh formula for Bessel functions. 



II. DARBOUX-BACKLUND TRANSFORMATIONS 



We begin to briefiy recall ther essential features of DBT and shape invariance. If ip\{x;a) is an eigenstate of the 
hamiltonian H{a) ~ —(P/dx^ + V{x;a), x e / C M, (a G R™ being a multiparameter) associated to the eigenvalue 
E),{a) {Eo{a) = 0) 
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^'i{x; a) + {Ex{a) - V{x; a)) ^x; a) = 0, (4) 

then the Riccati-Schrodinger (RS) function wx{x]a) = —ip')^{x;a)/ip\{x;a) satisfies the corresponding Riccati- 
Schrodinger (RS) equation^ 

- w'xix; a) + wl{x; a) = V{x; a) - Ex{a). (5) 

The class of RS equations is preserved by a specific subset of the group Q of smooth 5*^(2, ]R)-valued curves 
Map(R, SL{2,M.)) . These transformations, called Darboux-Backlund Transformations (DBT), are built from any 
solution Wu{x; a) of the initial RS equation Eq(l5]) as^"— 

wx(x;a) w\'(x;a)^-w„[x;a)-\ r r-, (6) 

w^(x;a) ~ wx[x;a) 

where Ex (a) ^ E^, (a) . w^'^^' is then a solution of the RS equation: 

- wi^^'ix; a) + (w[^'^ [x- a)) ' = F(^) (x; a) - Ex [a] , (7) 
with the same energy Ex{a) as in Eq(l5|) but with a modified potential 

V'-''\x;a) ^V{x;a) + 2wl{x;a) (8) 

called an extension of V. If V^'^^ is a rational function of an a priori given variable {x, tanx, sina;, coshx...) we call 
it a rational extension of V. 

This can be schematically summarized as 




(9) 



The corresponding eigenstate of H^'^^{a) — ~d? /dx^ + V'^^\x\a) can be written 

^ 0.P (- . ^=±==Mu,U.(.a, ,10, 

\ J ) ^Ex{a)-E^{a) 

where A (w^) is a first order operator given by A (wn) — d/dx + w^{x; a). Eq(ITU]) can still be written as 

where W {yi, ■■■,yn \ x) is the Wronskian of the functions j/i, 

From V , the DBT generates a new potential V'^"^ (quasi) isospectral to the original one and its eigenfunctions are 
directly obtained from those of V via Eq([TU|. Nevertheless, in general, 'Wy{x\a) and then the transformed potential 
is singular at the nodes of ij^^^x', a). For instance, if Eliji) admits a bound state spectrum (£'niV'n)„£N 

(x; a), F^"^ is regular only when n = 0, that is when ?/'n=o is the ground state of H, and we recover the usual SUSY 
partnership in quantum mechanicai^ii^. Note that A{wq) is a "state-deleting" DBT, the spectrum of the superpartner 

hamiltonian H = H^^^ having for fundamental level and ground state Ei and V'l"'' respectively. 
V is said to be a translationally shape invariant potential (TSIP) if its superpartner has the form 



V{x- a) = V^°\x- a) = V{x] a + a) + i?(a). 



(12) 
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where a G R"*. 

The question of successive iterations of DBT is very natural and is at the center of the construction of the hierarchy 
of hamiltonians in the usual SUSY QM scheme^. Staying at the formal level, it can be simply described by the 
following generalization of Eq® {Nj denotes the j-uple (y\, , i/j), with Ni = I'l and = {N,n-i, i^m)) 



W\ ^ W\ >— > ... w Wy^ 

where ui^^"^ is a RS function associated to the eigenvalue E\ of the potential 



(13) 



F^^") (x; a) = V{x- a) + 2 ^ (wif-^\x- a)) . (14) 



The corresponding eigenfunction is given by (cf Eq ([TU)) and Eg pT]) ') 



")(.;«) = A ^r'"-)(.;a) = ^ \, ' ' ^ ^ (15) 



or 



^^^-\x;a)^A(w\Jl--^) ...A{w.,)^j;x{x■,a). (16) 
By induction, we can prove the Crum formula a^'^'^^ 



V'l^'"''(a:;a) = (i/'^.i , , '0A | x) /W ....^^^ \ x) 
y(^'")(x; a) = V^(a:; a) - 2 (log W (^,, , | x))" . 



(17) 



A succession of m SUSY QM partnerships gives rise to a hierarchy of regular potentials j — l,...,m, 

associated to the m-uple Nm = (0, , m — 1) , that is, to a set of seeds functions {ipQ, ipm-i) constituted by the 
m first bound states of the initial hamiltonian. In the case where the initial potential is a TSIP, we obtain simply 

m — 1 

T/(^") {x- a) = V{x- a) - 2 (log W (V'o, -, V-m-i I x))" = V{x- a + ma) + ^ i?(a,), (18) 

with 



,,(Ar„)^^_^ _ W {lpo,...,i>m-l,'ipk+m \ x) 

W (ipo, ■.■,ipm-i I a;) 



^ZzHx^a)^ " r::: " ^ fc>o, as) 



which is the k-th excited state of y^^™) with the associated energy Ek+m- 

III. FROM THE CONSTANT POTENTIAL TO THE TDPT POTENTIAL 

Following Darboux^, our starting point is the constant potential on the real line V{x) = 0, x G K. The physical 
energy spectrum of the corresponding hamiltonian H is constituted by an open band £k = k"^ ^ ]0, +cxd[, fc > 0, 
where the two-dimensional eigenspace associated to £k is spanned by the (unnormalized) scattering eigenstates 

[x) = sin(fca;)/fc, 0^'"=" (x) = cos{kx)/k. (20) 
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We then build a chain of extensions of V where the first DBT we choose is based on a scattering state satisfying 
the Dirichlct boundary condition at the origin, namely 

^0 (x) = (pf'^ {x) = sin (x) , wq (x) = - cot (x) , (21) 
where the choice of the value A: = 1 is for pure convenience. The first extended potential is then 

V^°\x) = V{x) + 2w'o{x) = (22) 

sm X 

and the corresponding Schrodinger equation 

{^-cP/dx^ + y(°)(x)) i:{x) = Eil;{x). (23) 

At this step, Darbou»i writes 

"Pour des valeurs particulieres de E, cette equation admet les solutions sin^ x, sin^xcosx. En les employant et en 
poursuivant I 'application de la methode, on parviendra successivement d des equations de la forme 

,2/72 + 1) m(m + 1) \ , . s I / \ „ 

'dydx^ + „ H -ipix) = E'ipix)." 

COS'' X sin X / 

As emphasized by Gaillard and Matveev, the argument needs to be completed and specified in particular to 
determine in a systematic way what is the correct choice of the successive seeds functions to reach this equation. 

First note that Darboux's procedure does not consider the spectral properties of the built hamiltonians which can 
then become more singular at each step. This is the case of which is singular at every node = In of ipQ (x). 
The corresponding singularities are of the centrifugal barrier A/ (x — ^;)^ with A > 3/4, that is, strong singularities 
for which the transmission probability is zero^ii^ii^. This confines the particle in a single interval If we take 

for instance I = 0, the existence of the singularities in and tt can then be considered as equivalent to the addition 
of Dirichlet boundary conditions at the origin and at x = tt. 

ijj (0+) V (tt") = 0. (24) 

These lasts act as a filter which select among all the formal eigenf unctions ip^^ only those which satisfy Eq p4)l . 
The question is now to determine what are the appropriate values of k which allows to verify this constraint. 

From a more general point of view, consider the problem to determine what is the condition to which the solution 
4>\ (x) of EqQ is submitted in order that its image -01"^ via the DBT A (w^) satifies 

V'a(xo) = 0, (25) 
where xq € / is a zero of ip^. In the vinicity of xq we can write 

ip^,ix) = a^^o + a^a - ^o) + 0^,2 - ^of + O (^{x - xq)^^ . (26) 
Xq being a simple zero of ipi, we have a,y_o = and a^.i =/= 0, which gives 

w,{x) ~ -l]—-^ + 0{{x^xo)). (27) 

x-txo X ~ Xo Cti/,! 

Inserting these results in Eq ifTl]) . we obtain 

^M(^) ^ _ ^a,,o + 0((x - xo)) (28) 

x~>Xo X — Xo Ou 1 



and the condition (j25|) is satisfied iff a^^o = Oi that is, iff xo is also a zero of "0. 



A- 



5 



In brief, a DBT A {w^) associated to a seed function ijji, having nodes on the definition interval generates an 
extended hamihonian which is singular on this interval. The restriction of this extension between two singularities 
has a spectrum which consists in a filtered version of the initial hamiltonian spectrum: we only keep the levels E\ 
for which the initial eigenstates have nodes at the singular points. This point of view is equivalent to the one of 
Marquez et a£. 

Applying this analysis to to the system considered above, we obtain that the potential 1/^°' restricted to the positive 
half hne and subjected to the boundary conditions Eq ([24|) has a discrete non degenerate spectrum given by i?;, / G N*, 
with the corresponding eigenstates 

V'^"^ (x) = A {wq) iJi (x) ^ (l + l) cos ((; + 1) a;) - cot x sin {{I + I) x) , (29) 

where we have noted ipi (x) = ipf^lf (x) and Ei = £i+i- In particular, the fundamental and first excited eigenstates 

are ■i/'i"'' (x) = sin^ (x) and ipi^^ i^) = sin^ xcosx. V^'^'' can also be considered as the image via A (wq) of the infinite 
square well 



being a particular TDPT potential ([T]) with n = 0, m = 1, we can consider that Vsw correspond to the limit 
case Vsw{x) = V{x]0,0) = lim V{x; fi^v). 

As it is well knowiS^^—, the TDPT potentials V{x; ji, v) are second category translationally shape invariant 
potentials (TSIP), the SUSY QM partner of V{x] ^, v) being given by 



V{x; fi, I') — V{x; ^ + 1,1^+1), 



(31) 



for fj,,i> > and 



( V(x;fi,0) = V{x;ii + l,0), fi > 

\v{x;0,i^) = V{x;0,v + l), > 0. ^ ' 

Starting from V^'^^x) ~ V{x;l,0), we can build a hierarchy of p — 1 extensions of via successive SUSY QM 
partnerships, that is, via a chain of p — 1 DBT based on the successive ground states. From Eg ([32]) and Eq p8)) . we 
deduce that the final potential is V{x;p,0) and can be seen as the end of a chain of p extensions which starts from 
the zero potential V{x) and which is associated to the p-uple Np = {0, ...,p — 1), ie to the set of seeds eigenfunctions 

(-00 (x) , ipp-i {x)) = (sin (x) , sin {px)). Using a straightforward recurrence, we deduce that the ground state of 
yii'fp)(^x) = V{x;p, 0) = p{p +1)1 sin^ x at energy Ep is 

^fv) (2-) ^ sinP+i X. (33) 

Since is even with respect to 7r/2, its eigenfunctions ^i^"^ are repectively symmetric (if n—p is even) or antisymmetric 
(if n — p is odd) with respect to this point. An immediate consequence is that its first excited eigenstate (at energy 

Ep+i) V'p+i^ has only one node on ]0,7r[ which is always located at a; = 7r/2 independently of p. From Eg ([55]) we 
deduce 



4^11^^ {x) ^ A+ ^llt'^ {x) - cos X sinf +1 x. (34) 

Suppose that we continue the chain of extension and build the following DBT from this first excited eigenstate. 
The obtained potential is then singular and presents a pole of second order at a; = tt /2. Vp+i^ being the image of 
■0^+1(2;) = sin((p + 2) a;) by the chain associated to Np 



7f-^))...l(7«o)^P+i(a:), 



(35) 
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the new chain is associated to the (n + l)-uple {Np,p +1) and the final extension is given by 

= ^ + (P+^)^P + ^\ (36) 



2 

sm X 



The presence of the additional "strong" singularity in tt/2 imposes a Dirichlet boundary condition at this point. In 

the spectrum of the restriction of T/(^p'P+i) to ]0,7r/2[, half of the levels E^'''' of the preceding potential are 
eliminated. These are all those associated to "even" (with respect to 7r/2) eigenstates, that is to the even values of 
the quantum number n — p. 

The spectrum of y(^p'P+i) (x) contains then only the levels Ep^ij^2{j+i)i 3 ^ 0, and the associated eigenstates 
"0^+1+2(^+1) ^'^^ images of the initial eigenstates 0p+i+2(j+i) via the chain of BDT corresponding to the (p + 1)- 
uple {Np,p + 1) = (0, ...,p — l,p + 1). Starting from y(^p'P+i) we can continue the chain by using standard SUSY 
QM partnerships, all the extended potentials thus generated being perfectly regular on ]0, 7r/2[. But 1/(^p'P+i) is also 
a TDPT 

V^^^-P+'Hx) = V{x;p+l,l). (37) 
Consequently the shape invariance property Eq (|31|) implies that after q steps, we obtain as final extension 



V^'^'"P+^'-'P+i){x) = Vix;p + q,q). (38) 

This shows in an explicite way what is exactly the chain of DBT which permit to build the general TDPT potential 
of integer parameters as a rational (in sinx) extension of the free particle system. 
The Crum formula gives then 

V{x;p + q,q) = -2 (In (t/iq, ^p+i,---- V'p-i+2j, -0^-1+29 I x))" , (39) 

or, noting I — p + q, 

v( 1 \ g(g+l) , ^(^ + 1) 
cos^ X sm X 

= -2 (In W (sinx, sin((/ - q) x), sin ((/-<? + 2)a;) , sin ((/-(? + 2j)x) , sin ((/ + (?) a;) | x))" ,(40) 
which is precisely the Gaillard-Matveev result. 



IV. FROM THE CONSTANT POTENTIAL TO THE BESSEL POTENTIAL AND A WRONSKIAN 

RAYLEIGH FORMULA. 

To generate the TDPT potentials from the constant one, we start with a DBT built from any (physical) diffusion 
state of the constant potential. Nevertheless, we are not limited to this case and on a formal point of view we 
can use any eigenfunction, physical or not to build this DBT. In particular, we can choose as seed eigenfunction 
■i/'o (ic) = X = lim-i/ife [x) whose associated eigenvalue i?o = is located at the lower boundary of the physical spectrum 

of V . The extended potential generated by the DBT A{w^ is 



y(o)(x)-2/x2 = VB(a;,l), (41) 

that is, the first Bessel potential of integer parameter, the general form of which being Vb (x\ a) — a{a + 1) /x"^ . 

The two partner hamiltonians H and iJ^^-'are strictly isospectral but the levels of are no more degenerated since 
the strong sing ularity of at the origin imposes to retain in the spectrum of H only the eigenstates which satisfy 
the Dirichlet boundary condition at the origin 



i) (0+) = 0. 



(42) 
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The physical spectrum of i?^"' is then 



e'^^'^ = P e ]0,+oo[, fc > 

(^) — Ek-Eo ^(^ojV'fc (2;) = (— sin(fca;) + kxcos{kx)) /k^x, 



(43) 



the unphysical "fundamental" eigenfunction for fc = being ip^'^ [x) = ~ Vaml))^' {x). 



,,(0) 



We can start to build a chain of isospectral potentials by successive DBT based on the " fundamental" unphysical 
eigenfunctions of zero energy which are recessive at the origin (satisfying the Dirichlet boundary conditions at the 
origin). An immediate recurrence gives, for = (0, ...,0) 



m times 



(44) 



All the extensions V'^'^"^^ subject to the Dirichlet boundary condition at the origin are then strictly isospectral to 
with a unique energy band ]0, +oo[. These results are naturally well-known and this constitutes the first example 
of application that Darboux gave of his method in his first article on the subject-". They are also one of the key 
ingredient given by Duistermaat and Griinbaum^- in the proof that every trivial monodromy potential decaying at 
infinity can be considered as a rational extension of the free particle potential. 

Nevertheless, the question of the Wronskian representation of the extended potentials Ea Pl)) and of their eigen- 
functions is less trivial. Indeed, we are now in a case of confluency where all the successive DBT are built from 
eigenfunctions associated to the same energy. Then we cannot use the Crum formulas to express and its eigen- 

states in terms of Wronskians of distinct "seeds functions" (note that Aivj"^"^'^) is an annihilator for V'o^"'' 



and that t/'q^"^^'' cannot be considered as the image of V'o by the DBT A{w); '''^')). In this case, we can however 
use Matveev's formulas^Mi and express F^^™) and V'fc^'"'' in terms of " generalized Wronskians": 

Matveev's formulas: In the confluent case where the repeated DBT are built on eigenfunctions associated to the 
same value fcg of spectral parameter k {N„i = (fco, fcp) ), we can write 



m times 



V(^-Hx) ^V{x)-2 (\nW (^V., ^ 



K—ko 



■.Tpk\x 



(45) 
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where the 



K—ko 



j — 0, ...,m — 1, are the m first non-zero derivative of V'k with respect to the spectral 



parameter n at the value k = feg. In our case, for fcg — 0, we have 



^k (a:) = ^(-l)" 



n=0 



{2n+l): 



(46) 



The sequence of the ij is given by ij = 2j, j > with 



and the Matveev formula for the potential gives 



Vix;m) = 2 {\nW {x,x\ ..^x^""^ \ x))" . (48) 
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Darboux's result is readily recovered by using standard properties of the Wronskians-^ which allow us to write 

m(m-l)/2 



W{x,x^...,x'--' I x) =2;™ fgj Ty(l,y,...,y™-M2/)|^^,. =x™("+i)/22"(™-i)/2 (49) 

that is, y(^")(a;) = VBix;m) = TO(m + 1) /a;2. 

More interesting is the case of the eigenfunctions. Since "i/'i^"'' ^ solution of 

^"(x)+(fc2-^^^^^^)^(x) = 0, (50) 
satisfying the Dirichlet boundary condition at the origin, it can then be written as 

V'i^"' (x) ~ a;^/^J„+i/2(fcx) = y^xj,„(fca;), (51) 

where Jm+1/2 (a:) and jm {x) are the usual Bessel and spherical Bessel functions— 1^. 
But from the second Matveev's formula we also have 



,,(Ar,„) _ W {x,x^,...,x^"' \sm{kx)/k I x) 
M^(x,x3,...,a;2™-i | x) 



This can be considered as a Wronskian version of the Rayleigh formula for Bessel functions^. Indeed, the same 
handling as before^^ applied to the Wronskian at the numerator of Eg ([5^ gives 

W {x,x\...,X^"'-\sm{kx)/k I x) =2™(™+l)/2^(™+l)(™+2)/2^^^^y_ _^^ym-l^^(^-) I ^ ^53) 

where y — x^ and g{y) = sinc(fca;) = sin(fca;)/fcx. Then 



W {x,x^...,x^"'-\smikx)/k I x) = 2™("+i)/2^(™+i)("+2)/2 | gM(y) (54) 

and with Eq([49|. we finally obtain 

m+l / 1 rl \ ^ 

V.r'"^(.) = ^(^--j smc(fcx). (55) 
V. CONCLUSION 



In this letter we have shown how to recover in a simple way the content of the Gaillard-Matveev theorem which 
provides a Wronskian representation for the TDPT potentials. This is directly achieved by combining the Crum 
formula and shape invariance arguments applied to specific singular extensions of the constant potential. The same 
reasoning can be adapted to obtain Wronskian representations for the eigenfunctions of the Bessel potentials. This 
confluent case necessitates to employ Matveev's generalized Wronskians, obtaining then a Wronskian version of the 
Rayleigh formula. 

In2., Gaillard and Matveev also consider the case of the discrete Darboux-Poschl- Teller potentials (DDPT), for 
which they give Casoratian representation formulas. A treatment of these discrete systems by the preceding approach 
will be the object of further investigations. 
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